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Abstract 

A Hamiltonian formalism is used to describe ensembles of fields in terms of two canonically conju- 
gate functionals (one being the field probability density). The postulate that a classical ensemble is 
subject to nonclassical fluctuations of the field momentum density, of a strength determined solely 
by the field uncertainty, is shown to lead to a unique modification of the ensemble Hamiltonian. 
The modified equations of motion are equivalent to the quantum equations for a bosonic field, 
and thus this exact uncertainty principle provides a new approach to deriving and interpreting the 
properties of quantum ensembles. The examples of electromagnetic and gravitational fields are 
discussed. In the latter case the exact uncertainty approach specifies a unique operator ordering 
for the Wheeler-DeWitt and Ashtekar-Wheeler-DeWitt equations. 
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1. INTRODUCTION 



The Heisenberg uncertainty relation 



(1) 



for the rms position and momentum uncertainties of a quantum particle is a well known 
feature of quantum mechanics. It has recently been shown, however, that this relation is 
a consequence of a more fundamental connection between the statistics of complementary 
quantum observables. 

In particular, the distinguishing "nonclassical" property of complementary observables 
is that they cannot be simultaneously measured to an arbitrary accuracy. It is therefore 
natural to consider the decomposition of one such observable, the momentum say, into the 
sum of a "classical" and a "nonclassical" component: 



where the classical component, is defined as that observable closest to p (in a statistical 
sense) which is simultaneously measurable with the complementary observable x [1]. 

It turns out that such decompositions do indeed, in a number of ways, neatly separate clas- 
sical and nonclassical properties of quantum observables. For example, for one-dimensional 
particles the nonclassical component of the momentum satisfies the uncertainty relation 



for a// pure states, where 6x denotes a measure of position uncertainty from classical statistics 
called the Fisher length [1] . This exact uncertainty relation is far stronger than (and implies) 
the corresponding Heisenberg uncertainty relation in Eq. (1). 

The surprising fact that quantum particles satisfy exact uncertainty relations has recently 
provided the basis for deriving much of the quantum formalism, from an exact uncertainty 
principle. In particular, the assumption that a classical ensemble is subjected to nonclassical 
momentum fluctuations, of a strength inversely proportional to uncertainty in position, has 
been shown to lead directly from the classical equations of motion to the Schrodinger equa- 
tion [2]. A brief overview of exact uncertainty properties of quantum particles is provided 
by the conference paper in Ref. [3]. 



■p = p + p' 
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The aim of this paper is to show that the exact uncertainty principle may be successfully 
generalised to derive the equations of motion for bosonic fields with Hamiltonians quadratic 
in the field momenta (eg, scalar, electromagnetic and gravitational fields). This "exact un- 
certainty" approach is extremely minimalist in nature: unlike canonical quantisation, it does 
not use nor make any assumptions about the existence of operators, Hilbert spaces, complex 
amplitudes, inner products, linearity, superposition, or the like. The sole "nonclassical" 
element needed is the addition of fluctuations to the momentum density of members of a 
classical ensemble of fields, with the fluctuation statistics assumed to be determined by the 
ensemble statistics. The exact uncertainty approach is thus conceptually very simple, being 
based on the core notion of statistical uncertainty (intrinsic to any interpretation of quan- 
tum theory). As a bonus, the exact uncertainty approach further implies a unique operator 
ordering for the Schrodinger equation associated with the quantum ensemble - something 
which the canonical quantisation procedure is unable to do. 

It is remarkable that the basic underlying concept - the addition of "nonclassical" mo- 
mentum fluctuations to a classical ensemble - carries through from quantum particles to 
quantum fields, without creating conceptual difficulties (although significant technical gen- 
eralisations are needed). This logical consistency and range of applicability is a further 
strength of the exact uncertainty approach. 

In the next section the equations of motion for a classical ensemble of fields are expressed 
in Hamiltonian form, via two canonically conjugate functionals (the probability density 
and the Hamilton- Jacobi functional). In Sec. 3 it is shown that the classical ensemble 
Hamiltonian is modified by the addition of nonclassical momentum density fluctuations, in 
a manner uniquely specified by the exact uncertainty principle, leading to modified equations 
of motion equivalent to the quantum field equations. 

In Sees. 4 and 5 the examples of the electromagnetic and gravitational fields are dis- 
cussed. In the former case it is shown that the exact uncertainty approach is equivalent 
to adding nonclassical fluctuations to the electric held E, of a strength determined by the 
inherent uncertainty of the vector potential A. In the latter case it is shown that the ex- 
act uncertainty approach leads to a unique operator ordering for the Wheeler-DeWitt and 
Ashtekar-Whccler-DeWitt equations. This ordering is, moreover, consistent with Vilenkin's 
"tunneling" boundary condition for inflationary cosmology [4]. Results are discussed in 
Sec. 6, and necessary elements from classical fleld theory and functional analysis, and the 
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proof of the main Theorem of Sec. 3, are given in the Appendices. 



2. CLASSICAL ENSEMBLES 

We first consider a real multicomponent classical field / = (/") with conjugate momentum 
density g = {g""), described by some Hamiltonian functional H[f,g,t]. For example, / may 
denote the electromagnetic field A = (A^) , or some collection of interacting fields labelled 
by the index a. Spatial coordinates will be denoted by x (irrespective of dimension), and the 
values of field components and at position x will be denoted by and g^ respectively. 

There are three canonical approaches to describing the evolution of an individual field, 
based on the (related) Lagrangian, Hamiltonian, and Hamilton- Jacobi formalisms respec- 
tively. We choose the latter here, as it provides a straightforward mechanism for adding 
momentum fluctuations to an ensemble of fields (Eq. (8) below). The Hamilton- Jacobi for- 
malism also leads to an elegant "ensemble Hamiltonian" representation for the dynamics of 
an ensemble of fields, in terms of two canonically conjugate functionals (Eq. (5) below). 

First, the equation of motion for an individual classical field is given by the Hamilton- 
Jacobi equation 

— + H[f,6S/6f,t] = 0, (3) 

where S[f] denotes the Hamilton- Jacobi functional, and S/Sf denotes the functional deriva- 
tive with respect to / (see Appendices A and B). The momentum density associated with 
field / is given by = SS/Sf, and hence S will also be referred to as the momentum potential. 

Second, the description of an ensemble of such fields further requires a probability den- 
sity functional, P[f]- The equation of motion for P[f] corresponds to the conservation of 
probability, i.e., to the continuity equation 




as is reviewed in Appendix B. 

Eqs. (3) and (4) describe the motion of the ensemble completely, in terms of the two 
functionals P and S. These equations of motion may be written in the "Hamiltonian" form 



where H denotes the functional integral 

H[P, S, t] := {H) = jDf PH[f, 6S/6f, t], (6) 

and where variational derivatives such as AH/ AS are discussed in Appendix A. The equiv- 
alence of Eqs. (3) and (4) to Eqs. (5) follows directly from Eq. (A6). 

The functional integral H in Eq. (6) will therefore be referred to as the ensemble Hamil- 
tonian, and, in analogy to Eqs. (Bl) of Appendix B, P and S may be regarded as canonically 
conjugate functionals. Note from Eq. (6) that H typically corresponds to the mean energy 
of the ensemble. It may be shown that Eqs. (5) follow from the action principle AA — 0, 
with action A = / dt[-H + J Df S{dP/dt)]. 

In what follows, we specialise to ensembles for which the associated Hamiltonian func- 
tional is quadratic in the momentum field density, i.e., of the form 

H[f. g,t] = J2ldx KfiMgl + V[f] . (7) 

a,b 

Here K^'^ [f] = K^^ [f] is a kinetic factor coupling components of the momentum density, and 
V[f] is some potential energy functional. The corresponding ensemble Hamiltonian is given 
by Eq. (6). Note that cross terms of the form g^gt' with x ^ x' are not permitted in local 
field theories, and hence are not considered here. 

3. MOMENTUM FLUCTUATIONS ^ QUANTUM ENSEMBLES 

The ensemble Hamiltonian, H[P,S,t\ in Eq. (6), is our classical starting point for de- 
scribing an ensemble of fields. This starting point must be modified in some way if one is to 
obtain new equations of motion, to be identified as describing a quantum ensemble of fields. 
Our approach to modifying the ensemble Hamiltonian is based on a single ingredient: the 
addition of nonclassical fiuctuations to the momentum density, with the magnitude of the 
fiuctuations determined by the uncertainty in the field. This "exact uncertainty" approach 
leads to equations of motion equivalent to those of a bosonic field, with the interpretational 
advantage of an intuitive statistical picture for quantum field ensembles, and the technical 
advantage of a unique operator ordering for the associated Schrodinger equation. 

Suppose then that 6S/6f is in fact an average momentum density associated with field 
/, in the sense that the true momentum density is given by 

g^ SS/Sf + N, (8) 
5 



where is a fluctuation field that vanishes on the average for any given field /. Thus 
the physical meaning of S changes to being an average momentum potential. No specific 
underlying model for is assumed or necessary: in the approach to be followed, one may in 
fact interpret the "source" of the fluctuations as the fleld uncertainty itself. Thus the main 
effect of the fluctuation field is to remove any deterministic connection between / and g. 

Since the momentum fiuctuations may conceivably depend on the field /, the average 
over such fiuctuations for a given quantity A[f, N] will be denoted by A[f], and the average 
over fiuctuations anc? the field by (A). Thus = by assumption, and in general (A) — 
J Df P[f] A[f]. Assuming a quadratic dependence on momentum density as per Eq. (7), it 
follows that when the fluctuations are signiflcant the classical ensemble Hamiltonian H = 
{H) in Eq. (6) should be replaced by 

H' = {H[f,SS/Sf + N,t] ) 

a,b 

= H + J2jDf J dxPKfN^^. (9) 

a,b 

Thus the momentum fluctuations lead to an additional nonclassical term in the ensemble 
Hamiltonian, specifled by the covariance matrix Coyx{N) of the fluctuations at position x, 
where 

[Coy,(N)f -.^NSNI. (10) 

The covariance matrix is uniquely determined, up to a multiplicative constant, by the 
following four assumptions: 

(1) Causality: H' is an ensemble Hamiltonian for the canonically conjugate functionals P 
and which yields causal equations of motion. Thus no higher than first-order functional 
derivatives can appear in the additional term in Eq. (9), implying that 

Cov,(Ar) = a(P, 5P/5U, S, SS/Sf,, t) 

for some symmetric matrix function a. Note that in principle one could also allow the covari- 
ance matrix to depend on auxilary fields and functionals; however, the fourth assumption 
below immediately removes such a possibility. 

(2) Independence: If the ensemble comprises two independent non-interacting subensem- 
bles 1 and 2, with a factorisable probability density functional P[f^^\ f^^^ = Pi[f^^^]P2[f^% 
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then any dependence of the corresponding subensemble fluctuations N^^^ and N^"^^ on P only 
enters via the corresponding probability densities Pi and P2 respectively. Thus 

Cov.(7VW)|^^^^ = Cov.(A^^^^)|^^ , Cov.(7V(2))|^^^^ = Cov.(A^^^))|^^ 

for such an ensemble. Note that this assumption implies that the ensemble Hamiltonian 
H' in Eq. (9) is additive for independent non-interacting ensembles (as is the corresponding 
action A'). 

(3) Invariance: The covariance matrix transforms correctly under linear canonical trans- 
formations of the field components. Thus, noting that / — ^ A^^/, g A^g is a canonical 
transformation for any invertible matrix A with transpose A"^, which preserves the quadratic 
form of H in Eq. (7) and leaves the momentum potential S invariant (since 5/5 f ^ A^5 /5f), 
one has from Eq. (8) that N — > A'^N, and hence that 

CoY^{N) A^CoY^{N)A for / ^ A'V. 

Note that for single-component fields this reduces to a scaling relation for the variance of 
the fluctuations at each point x. 

(4) Exact uncertainty: The uncertainty of the momentum density fluctuations at any 
given position and time, as characterised by the covariance matrix of the fluctuations, is 
specifled by the fleld uncertainty at that position and time. Thus, since the fleld uncertainty 
is completely determined by the probabihty density functional P, it follows that CoVa;(A^) 
cannot depend on S, nor explicitly on t. 

It is seen that the flrst three assumptions (causality, independence and invariance) are 
natural on physical grounds, and hence relatively unconstraining. In contrast, the fourth 
assumption is "special" : it postulates an exact connection between the nonclassical momen- 
tum uncertainty and the fleld uncertainty. Remarkably, these assumptions lead directly to 
equations of motion of a bosonic quantum fleld, as shown by the following Theorem and 
Corollary. 

Theorem: The above assumptions of causality, independence, invariance, and exact 
uncertainty imply that 

NSNI = C{5P/5r^){5P/5fl)/P\ (11) 
where C is a positive universal constant. 
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The theorem thus yields a unique form for the additional term in Eq. (9), up to a 
multiplicative constant C. The classical equations of motion for the ensemble are recovered 
in the limit of small fluctuations, i.e., in the limit C ^ 0. Note that one cannot make the 
identification iV" ~ {5P/6f^)/P from Eq. (11), as this is inconsistent with the fundamental 
property A^^ = 0. The proof of the theorem is given in Appendix C, and is substantially 
different from (and stronger than) proofs of an analogous theorem for quantum particles 
[2, 3] (the latter proofs rely heavily on a "scalar" assumption that does not carry over in a 
natural manner to fields). 

The main result of this section is the following corollary (proved in Appendix C): 
Corollary: The equations of motion corresponding to the ensemble Hamiltonian H' can 
be expressed as the single complex equation 

where one defines 

h := 2VC, * := \fPe'^l^. (13) 

Eq. (12) may be recognised as the Schrodinger equation for a quantum bosonic field 
[5, 6], and hence the goal of deriving this equation, via an exact uncertainty principle for 
nonclassical momentum fluctuations acting on a classical ensemble, has been achieved. Note 
that the exact uncertainty approach specifies a unique operator ordering, {S/Sf^)K^^{S/Sf^), 
for the functional derivative operators in Eq. (12). Thus there is no ambiguity in the 
ordering for cases where K^'' depends on the field /, in contrast to traditional approaches 
(eg, the Wheeler-DeWitt equation, discussed in Sec. 5 below) . The above results generalise 
straightforwardly to complex classical flclds. 

The ensemble of fields corresponding to ensemble Hamiltonian H' will be called the 
quantum ensemble corresponding to H. Note from Eqs. (11) and (13) that the role of 
Planck's constant is to fix the relative scale of the nonclassical fiuctuations. It is remarkable 
that the four assumptions of causality, independence, invariance and exact uncertainty lead 
to a linear operator equation. 

Finally, it may be remarked that the equations of motion of a classical ensemble may 
be subject to some imposed constraint(s) on P and S. For example, each member of an 
ensemble of electromagnetic fields may have the Lorentz gauge imposed (see Sec. 4 below) . 
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As a guiding principle, we will require that the corresponding quantum ensemble is subject 
to the same constraint (s) on P and S. This will ensure a meaningful classical-quantum 
correspondence for the results of field measurements. However, consistency of the quantum 

equations of motion with a given set of constraints is not guaranteed by the above Theorem 
and Corollary, and so must be checked independently for each case. 

4. EXAMPLE: ELECTROMAGNETIC FIELD 
A. Lorentz gauge 

The electromagnetic field is described, up to gauge invariance, by a 4-component field 
A = (A^). In the Lorentz gauge all physical fields satisfy d^A^^ = 0, and the classical 
equations of motion in vacuum are given by d'^d^A'^ — 0. These follow, for example, from 
the Hamiltonian [6] 



where ri^„ denotes the Minkowski tensor, n'^ denotes the conjugate momentum density, and 
V denotes the spatial derivative. Here Hqb corresponds to the gauge- breaking Lagrangian 

L = —(1/2) / dx A^'^ A^^y, and is seen to have the quadratic form of Eq. (7) (with K^^ = 



The exact uncertainty approach therefore immediately implies, via the Corollary of the 
previous section, that the evolution of a quantum ensemble of electromagnetic fields is de- 
scribed by the Schrodinger equation 



in agreement with the Gupta-Blculcr formalism [6]. 

Further, note that the probability of a member of the classical ensemble not satisfying 
the Lorentz gauge condition = is zero by assumption, i.e., the Lorentz gauge is 

equivalent to the condition that the product {dnA^)P[A] vanishes for all physical fields. For 
the quantum ensemble to satisfy this condition, as per the guiding principle discussed at the 
end of Sec. 3 above, one equivalently requires, noting Eq. (13), that 




(14) 



in{d^/dt) = Hgb[A, -in{5/5A)]m, 



(15) 



{d^A^mA]^Q. 
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As is well known, this constraint, if initially satisfied, is satisfied for all times [7] (as is 
the weaker constraint that only the 4-divergence of the positive frequency part of the field 
vanishes [6]). Hence the evolution of the quantum ensemble is consistent with the Lorentz 
gauge. It would be of interest to derive the consistency of this constraint directly from the 
equations of motion, Eqs. (CI) and (C2) of Appendix C, for P and S. 

B. Radiation gauge 

It is well known that one can also obtain the classical equations of motion for the electro- 
magnetic field via an alternative Hamiltonian, obtained by exploiting the degree of freedom 
left by the Lorentz gauge to remove a dynamical coordinate (corresponding to the longitu- 
dinal polarisation). In particular, since d^A^ is invariant under ^ + d^x for any 
function x satisfying d'^d^x = Oj one may completely fix the gauge in a given Lorentz frame 
by choosing x such that AP — 0. One thus obtains, writing = (A°,A), the radiation 
gauge A° = 0, V • A = 0. 

The classical equations of motion, d'^d^.A — 0, follow from the Hamiltonian 

i/fl[A,7r] = (l/2) ydx(7r-7r/eo+|Vx A|V//o), (16) 

where tt denotes the conjugate momentum density Here Hr corresponds to the standard 
Lagrangian L = j dxP^'^F^^, with A° = 0. 

This Hamiltonian has the quadratic form of Eq. (7), and hence the exact uncertainty 
approach yields the corresponding Schrodinger equation 

ih{d^/dt) = Hr[A, -in{5/5A)\^ (17) 

for a quantum ensemble of electromagnetic fields in the radiation gauge (this is also the form 
of the Schrodinger equation obtained via the Schwinger-Tomonaga formahsm [8]). 
Note that the electric field follows via Eqs. (16) and (Bl) as 

E = -dA/dt = -SHr/Stt = -n/eo, 

and is therefore directly proportional to the classical momentum density tt. Fluctuations of 
the momentum density thus correspond to fiuctuations of the electric field E. Further, the 
constraint V • A = implies there is a one-one relation between A and the magnetic field 
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B = V X A (up to an additive constant) . Uncertainty in the vector potential thus corresponds 
to uncertainty in the magnetic field B. Hence, in the radiation gauge, the exact uncertainty 
approach corresponds to adding nonclassical fiuctuations to the electric field components 
of an ensemble of electromagnetic fields, with the fluctuation strength determined by the 
uncertainty in the magnetic field components. 

5. EXAMPLE: GRAVITATIONAL FIELD 
A. Hamilton-Jacobi constraints 

The gravitational field is described, up to arbitrary coordinate transformations, by the 
metric tensor g = {Qhv)- The corresponding invariant length may be decomposed as [9] 

ds^ = Qi^^dx^dx" = -(a^ - /3 • f3)dt'^ + 2Pidx'dt + -fijdx'dx^, 

in terms of the lapse function a, the shift function (3, and the spatial 3-metric 7 = (7^^). 
The equations of motion are the Einstein field equations, which follow from the Hamiltonian 
functional [9] 

H[y,7i,a,(3] = j dxaHGi'y,^^] "2 j dx (3iTi'^ ^. , (18) 

where tt = {'k^-') denotes the momentum density conjugate to 7, \j denotes the covariant 
3-derivative, and the Hamiltonian density Hq is given by 

Hg = (l/2)G,,-fc^[7]7r^%^' - 2 (%[7](det7)'/'. (19) 

Here ^^^R is the curvature scalar corresponding to 7ij, and 

Gijkib] = {liklji + liiljk - 7ii7w)(det7)~^/^. 

The Hamiltonian functional H corresponds to the standard Lagrangian L — 
J dx {— det gY^^Rlg], where the momenta vr*^ and tt* conjugate to a and Pi respectively 
vanish identically. However, the lack of dependence of H on 7r° and tt* is consistently main- 
tained only if the rates of change of these momenta also vanish, i.e., noting Eq. (Bl) of 
Appendix B, only if the constraints [9] 

5H/Sa = Tic = 0, SH/5(3i = -27r'^| • = (20) 
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are satisfied. Thus the dynamics of the field is independent of a and /3, so that these 
functions may be fixed arbitrarily. Moreover, these constraints immediately yield if = in 
Eq. (18), and hence the system is static, with no explicit time dependence. 

It follows that, in the Hamilton- J acohi formulation of the equations of motion (see Ap- 
pendix B), the momentum potential S is independent of a, /3 and i. Noting that tt = 5S/5^ 
in this formulation, Eqs. (20) therefore yield the corresponding constraints 

_ 5S_ _ dS_ _ (6S_\ 

for S. As shown by Peres [10], a given functional FYf\ of the 3-metric is invariant under 
spatial coordinate transformations if and only if {6F/S^ij)\j = 0, and hence the fourth 
constraint in Eq. (21) is equivalent to the invariance of S under such transformations. This 
fourth constraint moreover implies that the second term in Eq. (18) may be dropped from 
the Hamiltonian, yielding the reduced Hamiltonian 

Ha['-f,TT,a] — j dxaHG["y,T^] (22) 

in the Hamiltonian- Jacobi formulation [10, 11]. 

For an ensemble of classical gravitational fields, the independence of the dynamics with 
respect to a, /3 and t implies that members of the ensemble are distinguishable only by 
their corresponding 3-metric 7. Moreover, it is natural to impose the additional geometric 
requirement that the ensemble is invariant under spatial coordinate transformations. One 
therefore has the constraints 

6P 6P dP ^ f 6P\ 

for the corresponding probabihty density functional P[^\, analogous to Eq. (21). The first 
two constraints imply that ensemble averages only involve integration over 7. 

B. Quantum ensembles and operator-ordering 

Noting Eq. (19), the Hamiltonian Hq in Eq. (22) has the quadratic form of Eq. (7). Hence 
the exact uncertainty approach is applicable, and immediately leads to the Schrodinger 
equation 

ihd-^/dt = j dxaUab, -iH5/5^)]^ (24) 
12 



for a quantum ensemble of gravitational fields, as per the Corollary of Sec. 3. 

As discussed at the end of Sec. 3, we follow the guiding principle that all constraints 
imposed on the classical ensemble should be carried over to corresponding constraints on 
the quantum ensemble. Thus, from Eqs. (21) and (23) we require that P and S*, and 
hence in Eq. (13), are independent of a, /3 and t and invariant under spatial coordinate 
transformations, i.e.. 

Applying the first and third of these constraints to Eq. (24) immediately yields, via Eq. (19), 
the reduced Schrodinger equation 

Hob, -in{5/5^)]^ = (-/,72)^G,,,,[7]^* - 2 (=^)i?[7](det7)'/'* = 0, (26) 

which may be recognised as the Wheeler- De Witt equation in the metric representation [9] . 

A notable feature of Eq. (26) is that the Wheeler- De Witt equation has not only been 
derived from an exact uncertainty principle: it has, as a consequence of Eq. (12), been 
derived with a precisely defined operator ordering (with Gijki sandwiched between the two 
functional derivatives). Thus the exact uncertainty approach does not admit ambiguity in 
this respect, unlike the standard approach [9]. Such removal of ambiguity is essential to 
making definite physical predictions, and hence may be regarded as an advantage of the 
exact uncertainty approach. 

For example, Kontoleon and Wiltshire [12] have pointed out that Vilenkin's prediction of 
infiation in minisupcrspacc, from a corresponding Wheeler-DeWitt equation with "tunnel- 
ing" boundary conditions [4], depends critically upon the operator ordering used. In par- 
ticular, considering the class of orderings defined by an integer power p, with corresponding 
Wheeler-DeWitt equation [4] 
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* = (27) 



da^ a da a^ dcfP' 

(for a Fricdmann- Robertson- Walker metric coupled to a scalar field (p), Kontoleon and Wilt- 
shire show that Vilenkin's approach fails for orderings with p > \ [12]. Moreover, they 
suggest that the only natural ordering is in fact the "Laplacian" ordering corresponding to 
p— 1-1 which has been justified on geometric grounds by Hawking and Page [13]. 

However, noting that the relevant Hamiltonian functional in Eq. (2.7) of Ref. [4] is 
quadratic in the momentum densities of the metric and the scalar field, the exact uncertainty 
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approach may be applied, and yields the Wheeler-DeWitt equation corresponding to p = —1 
in Eq. (27). Hence the criticism in Ref. [12] is avoided. One also has the nice feature that 
the associated Wheeler-DeWitt equation can be exactly solved for this "exact uncertainty" 
ordering [4]. 

A certain degree of ambiguity remains, which derives from the need to introduce some sort 
of regularisation scheme to remove divergences arising from the product of two functional 
derivatives acting at the same point in the Wheeler-DeWitt equation. Such considerations, 
however, do not play a role in the example that we have just discussed, which concerns min- 
isuperspace quantisation involving a finite number of degrees of freedom. It is important to 
distinguish this regularisation problem from the far more difficult one associated with the 
requirement of "Dirac consistency", i.e., the need to find a choice of operator ordering and 
regularisation scheme that will permit mapping the classical Poisson bracket algebra of con- 
straints to an algebra of operators within the context of the Dirac quantisation of canonical 
gravity [14] . Our approach is based on the Hamilton- Jacobi formulation of classical gravity 
and, as shown by Bergmann [15], the functional form of the Hamilton- Jacobi functional ^S" 
is already invariant under the action of the group generated by the constraints. 

Finally, note that a similar approach may be applied to the Ashtckar formalism for gravi- 
tational fields [16] , where again the Hamiltonian is quadratic in the field momentum density 
(in particular, the two constraints linear in the momentum density a become constraints on 
the Hamilton- Jacobi functional S, corresponding to the invariance of S under spatial coor- 
dinate and internal gauge transformations [17] , while the constraint quadratic in a generates 
the Ashtekar- Wheeler-DeWitt equation with a unique operator ordering). 

6. DISCUSSION 

The main result of this paper is the derivation of the quantum equation of motion, Eq. 
(12), from an exact uncertainty principle, for fields with Hamiltonian functionals quadratic 
in the momentum density. 

It is important to emphasise that the exact uncertainty approach does not assume the ex- 
istence of a complex amplitude functional ^[/], nor the representation of fields by operators, 
nor the existence of a universal constant h with units of action, nor the existence of a linear 
operator equation in some Hilbert space. Only the assumptions of causality, independence. 
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invariance and exact uncertainty are required, all formulated in terms of a single nonclassical 
element (the uncertainty introduced by the momentum fluctuation N) . Since uncertainty is 
at the conceptual core of quantum mechanics, this is an elegant and pleasing result. 

The assumptions used also provide an intuitive picture for the origin of the Schrodingcr 
equation for bosonic fields, as arising from nonclassical fluctuations of the momentum den- 
sity. Of course this picture has limitations - the fluctuations arise from the uncertainty 
of the field itself, rather than from some external source, and hence are most certainly 
"nonclassical" rather than "classical" in nature. 

A minimalist interpretation of the exact uncertainty approach, based on Eqs. (8) and 
(11), is as follows. Every physical field has an intrinsic uncertainty, which is modelled by a 
corresponding statistical ensemble. Further, the nature of this inherent uncertainty is such 
as to preclude a deterministic relationship between the field and its conjugate momentum 
density - one must introduce fluctuations into the classical relationship, as per Eq. (8). 
However, the degree of indeterminism in this relationship is precisely quantifiable, in a 
statistical sense, being directly specified by the inherent field uncertainty as per Eq. (11). 

The above interpretation may be regarded as a significant sharpening of the so-called 
"statistical interpretation" of quantum mechanics [18], and is notably very different to the 
"causal interpretation" of Bohm and co-workers [19]. In the latter it is assumed that there 
is a pre-existing complex amplitude functional ^[/] = y/P exp(iS/fi) obeying a Schrodinger 
equation, which acts upon a single classical field via the addition of a "quantum potential" , 
Q[P], to the classical Hamiltonian. It is further assumed that the momentum density is pre- 
cisely g = 5S/6f, and that physical ensembles of fields have probabihty density functional 
P = l^p. In contrast, the exact uncertainty approach does not postulate the existence of 
adjunct amplitudes and potentials; the Schrodinger equation directly represents the evolu- 
tion of an ensemble, rather than of an external amplitude functional acting on individual 
systems (and is derived rather than postulated); and the basic tenet in Eq. (8) is that 
g ^ 5S/5f. 

A further strength of the exact uncertainty approach is that the basic underlying con- 
cept - the addition of "nonclassical" momentum fluctuations to a classical ensemble - carries 
through from quantum particles to quantum fields without creating conceptual difficulties. 
This adds an interpretational strength to the exact uncertainty approach not mirrored in 
other approaches that rely on connecting the equations of motion of classical and quan- 
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turn ensembles. For example, the abovementioned "causal interpretation" of Bohm and 
co-workers is explicity non-local, and hence non-causal, for relativistic fields [19]. As an- 
other example, one cannot simultaneously describe both the electric and the magnetic fields 
in generalisations of Nelson's stochastic approach to electromagnetic fields [20] . 

It is of interest to consider the scope and limitations of the exact uncertainty approach 
to physical systems. This approach has previously been applied to quantum particles [2, 3], 
and may be generalised to obtain the Pauli equation for a non-relativistic spin-1/2 particle, 
and the Schrodinger equation for particles with position-dependent mass (where in the lat- 
ter case one obtains the unique ordering ■p[2m{x)]~^p + V{x) for the Hamiltonian operator, 
corresponding to the ordering parameter a = in Rcf. [21]). In this paper the approach 
has been further generalised to bosonic quantum fields with Hamiltonians quadratic in the 
momentum density (including all relativistic integer-spin fields). It is also, indirectly, ap- 
plicable to the non-quadratic Hamiltonian functional of a nonrelativistic Schrodinger field 
(corresponding to second quantisation of the particle Schrodinger equation), in the sense 
that this field may be obtained as a low-energy limit of the complex Klein-Gordon field [5] 
(to which the exact uncertainty approach directly applies). 

However, a major question to be addressed in the future is whether the exact uncertainty 
approach is applicable to the derivation of fermionic field equations. These have two features 
which present challenges: the corresponding ensemble Hamiltonian is usually linear in the 
momentum density, and the anticommutation relations make it difficult to connect the 
equations of motion with corresponding classical equations of motion in the limit as — > 0. 
One possible approach is to determine whether exact uncertainty relations exist for such 
fields, analogous to Eq. (2), as these might suggest the statistical properties required by 
suitable "nonclassical" fluctuations (note that the exact uncertainty relations satisfied by 
bosonic fields are derived in Ref. [22]). 

Finally, in this paper the basic Schrodinger equation for bosonic fields has been obtained, 
with the advantageous features of an intuitive picture for the origin of the "quantum" nature 
of such fields, and a unique operator ordering in cases where other approaches are ambiguous. 
It would be of interest to consider further issues, such as the representation of general 
physical obscrvables by operators (addressed for the case of particles in Ref. [2]), boundary 
conditions, infinities, etc, from the new perspective on the conceptual and technical basis of 
quantisation offered by the exact uncertainty approach. 
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APPENDIX A: FUNCTIONAL DERIVATIVES AND INTEGRALS 



The necessary definitions and properties of functionals are noted here, including varia- 
tional properties of functional integrals. 

A functional, F[f], is a mapping from a set of physical fields (assumed to form a vector 
space) to the real or complex numbers, and the functional derivative of F[f] is defined via 
the variation of F with respect to /, i.e., 

5F = F[f + 6f] - F[f] = Jdx^^ 5U (Al) 

for arbitrary infinitesimal variations 5f . Thus the functional derivative is a field density, 
SF/Sf, having the value 5F/Sfx at position x. For curved spaces one may explicitly include 
a volume element in the integral, thus redefining the functional derivative by a multiplicative 
function of x; however, this is merely a matter of taste and will not be adopted here. The 
functional derivative is assumed to always exist for the functionals in this paper. 

It follows directly from Eq. (Al) that the functional derivative satisfies product and 
chain rules analogous to ordinary differentiation. The choice F[f] — f^' in Eq. (Al) yields 
5 fx' fx = 5{.x — x'). Moreover, if the field depends on some parameter, t say, then choosing 
5f. = f.{t + 6t) - fS) in Eq. (Al) yields 

dF dF /• , 5F df^ , , , 

for the rate of change of F with respect to t. 

Functional integrals correspond to integration of functionals over the vector space of 
physical fields (or equivalence classes thereof). The only property we require for this paper 
is the existence of a measure Df on this vector space which is translation invariant^ i.e., 
/ Df = J Df for any translation f' — f + h (which follows immediately, for example, from 
the discretisation approach to functional integration [5]). In particular, this property implies 
the useful result 

jDf— = for J DfF[f] <oo, (A3) 

which is used repeatedly below and in the text. Eq. (A3) follows by noting that the finiteness 
condition and translation invariance imply 

O^Jof {F[f + Sf] - F[f]) ^ jdx5f,^j Df 6F/5f, 
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for arbitrary infinitesimal translations. 

Thus, for example, if F[f] has a finite expectation value with respect to some probability 
density functional -P[/], then Eq. (A3) yields the "integration by parts" formula 

I DfP(5F/5f) ^-j Df(5P/5f)F. 

Moreover, from Eq. (A3) the total probabihty, / Df P, is conserved for any probabihty fiow 
satisfying a continuity equation of the form 




providing that the average flow rate, {V^), is flnite. 
Finally, consider a functional integral of the form 

/[F] = j Dfi{F,5F/5f), (A5) 

where ^ denotes any function of some functional F and its functional derivative. Variation 
of I[F] with respect to F then gives, to flrst order. 



A/ = I[F + AF] - I[F] ^ jof |(ae/aF)AF j dx [di/d{5F/5U)] [S{AF)/Sf, 

^ jof ^{di/dF) - J dx-^m/di6F/6f,)]^ AF 
^jdxjDfj- {[di/d{5F/5f,)] AF} . 



Assuming that the functional integral of the expression in curly brackets in the last term is 
finite, this term vanishes from Eq. (A3), yielding the result 

A/ 



analogous to Eq. (Al), where the variational derivative A//AF is defined by 



- dF r''6L 



d{SF/Sf,) 



(A6) 



AF OF J 5U 

A similar result holds for multicomponent fields, with summation over the discrete index a 
in the second term. 
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APPENDIX B: HAMILTON-JACOBI ENSEMBLES 



The salient aspects of the Hamilton- Jacobi formulation of classical field theory [23] are 
collected here, with particular attention to the origin of the associated continuity equation 
for ensembles of classical fields, required in Sees. 2 and 3. 

Two classical fields /, g are canonically conjugate if there is a Hamiltonian functional 
H[f,g, t] such that 

df/dt = SH/Sg, dg/dt = -SH/6f. (Bl) 

These equations follow from the action principle 5A — 0, with action functional A — 
J dt [—H + / dx gx{dfx/dt)]. The rate of change of an arbitrary functional G[f,g,t\ follows 
from Eqs. (A2) and (Bl) as 

^ - ^ + [dx (^^JL _ ^^j£\ _. ^ + |(^ ^1 
dt dt J \5f^Sg^ 5g^5f^) ' dt 

where { , } is a generalised Poisson bracket. 

A canonical transformation maps /, g and H to /', g' and H' ^ such that the equations of 
motion for the latter retain the canonical form of Eq. (Bl). Equating the variations of the 
corresponding actions A and A' to zero, it follows that all physical trajectories must satisfy 

-H + j dx g:c{dfxldt) = -H' + J dx g'^{df!^/dt) + dF/dt 

for some "generating functional" F. Now, any two of the fields f,g,f',g' determine the 
remaining two fields for a given canonical transformation. Choosing / and g' as the two 
independent fields, defining the new generating functional G[f,g',t] — F + J dx f!^g'^, and 
using Eq. (A2), then yields 



+ / dx 



dt \5U ) dt \5g'^ 

for all physical trajectories. The terms in round brackets therefore vanish identically, yielding 
the generating relations 

H' ^ H + dG/dt, g = SG/Sf, f' = 6G/Sg'. (B2) 

A canonical transformation is thus completely specified by the associated generating func- 
tional G. 
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To obtain the Hamilton- Jacobi formulation of the equations of motion, consider a canon- 
ical transformation to fields g' which are time-independent (eg, to the fields / and g 
at some fixed time to). Prom Eq. (Bl) one may choose the corresponding Hamiltonian 
H' = without loss of generality, and hence from Eq. (B2) the momentum density and the 
associated generating functional S are specified by the functional equations 

^ = ^ + H[f,5S/Sf,t]=0. (B3) 

The latter is the desired Hamilton- Jacobi equation. Solving this equation for S is equivalent 
to solving Eqs. (Bl) for / and g. 

Note that along a physical trajectory one has g' = constant, and hence from Eqs. (A2) 
and (B3) that 

dS dS SS df^ „ df^ dA 



dt dt J Sf^ dt J dt dt 

Thus the Hamilton-Jacobi functional S is equal to the action functional A, up to an additive 
constant. This relation underlies the connection between the derivation of the Hamilton- 
Jacobi equation from a particular type of canonical transformation, as above, and the deriva- 
tion from a particular type of variation of the action, as per the Schwinger-Tomonaga for- 
mahsm [8, 24]. 

The Hamilton-Jacobi formulation has the interesting feature that once S is specified, the 
momentum density is determined by the relation g = dS/df, i.e., it is a functional of /. 
Thus, unlike the Hamiltonian formulation of Eqs. (Bl), an ensemble of fields is specified by 
a probability density functional P[f], not by a phase space density functional p[f,g]- 

In either case, the equation of motion for the probability density corresponds to the 
conservation of probability, i.e., to a continuity equation as per Eq. (A4). For example, in 
the Hamiltonian formulation the associated continuity equation for p[f, g] is 

dp/dt + I dx{{5/5U)[p{dU/dt)] + {5/5gMd9x/dt)]} = 0, 

which reduces to the Liouville equation dp/dt — {H,p} via Eqs. (Bl). 

Similarly, in the Hamilton-Jacobi formulation, the rate of change of the field / follows 
from Eqs. (Bl) and (B3) as the functional 



V,[f] = dfjdt = {5H/5g,) 



9=5S/Sf 
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and hence the associated continuity equation for an ensemble of fields described by P[f] 
follows via Eq. (A4) as 



Eqs. (B3) and (B4) generalise immediately to multicomponent fields, and form the basis 
of the classical starting point in the derivation of the quantum equations of motion for 
bosonic fields in Sees. 2 and 3. 

APPENDIX C: PROOFS OF THE THEOREM AND COROLLARY 

Proof of Theorem (Eq. 11): Prom the causality and exact uncertainty assumptions 
in Section 3, one has CoVa;(A'") = a{P,5P/5fx, fx)- To avoid issues of regularisation, it is 
convenient to consider a position-dependent canonical transformation, — > ^x^fx, such 
that A[A] := exp[/ dx In | det Aa;|] is finite. Then the probability density functional P and 
the measure Df transform as P ^ AP and Df — > A~^Df respectively, and so the invariance 
assumption in Section 3 requires that 



where u"' and denote the vectors SP/Sf^ and respectively, for a given value of x. 
Since A^. can remain the same at a given point x while varying elsewhere, this homogeneity 
condition must hold for A and A^ independently. Thus, choosing A^; to be the identity 
matrix at some point x, one has a{AP, Au,w) = a{P,u,w) for all A, implying that a can 
involve P only via the combination v :— u/P. 

The above homogeneity condition for a therefore reduces to 



Note that this equation is linear, and invariant under multiplication of a by any function of 

the scalar J := v'^w. Moreover, it may easily be checked that if a and r are solutions, then 
so are ar'^a and Ta~^T. Choosing the two independent solutions a = vv^, r = {ww^)~^, 
it follows that the general solution has the form 




(B4) 



a{AP, AA^u, A^ ^w) = A^q;(P, u, w)A, 



■xj 



1 



for arbitrary functions (3 and 7. 
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For P — P1P2 one finds v — {vi,V2).i w — (u^i, wq), where the subscripts label correspond- 
ing subensemble quantities, and hence the independence assumption in Section 3 reduces to 
the requirements 

/3(Ji + J2)vivl + 7(Ji + J2)KO-' = PiiJi)vivJ + 7i(Ji)(«;iO-\ 

P{Jl + J2)V2VI + 7(^1 + J2){W2WI)-^ = (32{J2)V2VJ + -f2{J2){w2W^)-\ 

for the respective subensemble covariance matrices. Thus /? = /3i = /32 = C*, 7 = 7i = 72 = 
D for universal (i.e., system- independent) constants C and D, yielding the general form 

[Cov.(Ar)]"^ = C{5P/5f:){5P/8fl)/P^ + DWf[f] 

for the fluctuation covariance matrix, where denotes the inverse of the matrix with 

a6-coefiicient /"/^. 

Since Wx [/] is purely a functional of /, it merely contributes a classical additive potential 
term to the ensemble Hamiltonian of Eq. (9). It thus has no nonclassical role, and can 
be absorbed directly into the classical potential iy) (indeed, for fields with more than 
one component this term is singularly ill-defined, and hence can be discarded on physical 
grounds). Thus we may take D — without loss of generahty. Finally, the positivity of C 
follows from the positivity of the covariance matrix Cov^(A'"), and the theorem is proved. 

Proof of Corollary (Eq. 12): First, the equations of motion corresponding to the 
ensemble Hamiltonian H' follow via the theorem and Eqs. (5) as: (a) the continuity equation 
Eq. (4) as before (since the additional term does not depend on S), which from Eq. (7) has 
the explicit form 

and (b) the modified Hamilton- Jacobi equation 

dS/dt = -AH'/AP = -H[f, 6S/6f, t] - A{H' - H)/AP. 
Calculating the last term via Eq. (11) and Eq. (A6) of Appendix A, this simphfies to 

Second, writing ^ = P^^'^ ex.p{iS/h), multiplying each side of Eq. (12) on the left by 
"if^, and expanding, gives a complex equation for P and S. The imaginary part is just the 
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continuity equation of Eq. (CI), and the real part is the modified Hamilton- Jacobi equation 
of Eq. (C2) above, providing that one identifies C with /A. 
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